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Abstract 

> 

^vq ■ We express the amplitudes for charmless three-body B decays in terms of diagrams. 

In addition, we show how to use Dalitz-plot analyses to obtain decay amplitudes 
^ . which are symmetric or antisymmetric under the exchange of two of the final-state 

particles. When annihilation-type diagrams are neglected, as in two-body decays, 
many of the exact, purely isosp in-based results are modified, leading to new tests of 
the standard model (SM). Some of the tests can be performed now, and we find that 
present data agree with the predictions of the SM. Furthermore, contrary to what 
was thought previously, it is possible to cleanly extract weak-phase information from 
three-body decays, and we discuss methods for B — y Ktctc, KKK, KKtt and ititit. 
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1 Introduction 



The B-factories BaBar and Belle ran for over ten years, and made an enormous 
number of measurements of observables in B decays. For the most part, these 
decays were of the form B — > MiM 2 (Mi is a meson), as these are most accessible 
experimentally. Nevertheless, there have still been some probes of three-body B — > 
MiM 2 M 3 decays. To be specific, experiments have obtained Dalitz plots for many 
of the decay modes in B — > Ktttt, KKK, KKtt, tttttt, and made measurements of 
(or obtained upper limits on) the branching ratios and indirect (mixing-induced) 
CP asymmetries of a number of these decays pQ. 

Things are similar on the theory side. The vast majority of theoretical analyses 
involve two-body B decays. This is in part due to the relative angular momentum 
of the final-state particles. For example, consider B® — > tt + tt~. Because there are 
two particles in the final state, it has a fixed value of I (in this case 1 = 0), and 
so tt + tt~ is a CP eigenstate. On the other hand, in the decay B% — > Kstt + tt~ , the 
7r + 7r~ can have even or odd relative angular momentum, so that Kstt + tt~ is not 
a CP eigenstate. This makes it much more difficult to find clean predictions of 
the standard model (SM) to compare with experimental measurements. This is a 
general property of three-body decays. 

Still, there have been some theoretical analyses of CP-conserving observables in 
three-body B — > Ktttt, KKK decays (2J [3j HJ |5] . In general, these studies examined 
the isospin decomposition of the decay amplitudes, and symmetry relations among 
them. The analyses were carried out using isospin amplitudes. 

In this paper, we examine the amplitudes of the three-body charmless decays 
B — > Ktttt, KKK, KKtt, tttttt using diagrams. In addition, using Dalitz-plot anal- 
yses of such decays, we show how to separate the amplitudes into pieces which 
are symmetric or antisymmetric under the exchange of two of the final-state parti- 
cles. This is useful for any decay which contains particles which are identical under 
isospin. Now, as has been shown in Ref. [Bj, the amplitudes for two-body B decays 
can be expressed in terms of 9 diagrams. However, 3 of these - the annihilation- 
type diagrams - are expected to be quite a bit smaller than the others, and can 
be neglected, to a good approximation. This same procedure can be applied to 
three-body decays. 

The point of this is as follows. When one neglects annihilation-type diagrams, 
new features appear. A given set of three-body decays (e.g. B — > Ktttt) contains a 
number of different transitions (e.g. B + — > K + tt + tt~, B\ — > K + tt°tt~, etc.). There 
are exact relations among the symmetric or antisymmetric amplitudes for these 
specific decays. However, when one neglects certain diagrams, these relations can 
be modified, and this can lead to new effects. For example, some linear combinations 
of the isospin amplitudes vanish for certain decays. Also, there are additional tests 
of the SM. In some cases, it is even possible to obtain clean information about the 
CP-violating phases. 
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In Sec. 2, we present the diagrams describing B — > M1M2M3 processes. We 
review Dalitz-plot analyses of three-body decays in Sec. 3, and show how to obtain 
amplitudes which are symmetric or antisymmetric under the exchange of two of the 
final-state particles. The decays B — > Kirn, B — > KKK, B — > KKtc and B — > -rhh 
are discussed in Sees. 4, 5, 6 and 7, respectively. In all cases, we give the expressions 
for the decay amplitudes in terms of diagrams, and examine the prospects for the 
clean extraction of weak-phase information. Other subjects related to the particular 
decays are also discussed: resonances and penguin dominance in B — > Kirn (Sec. 4), 
penguin dominance and isospin amplitudes in B — > KKK (Sec. 5), T dominance in 
B — > KKix (Sec. 6), and Dalitz plots in B — > -nun (Sec. 7). We conclude in Sec. 8. 

2 Diagrams 

It has been shown in Ref. [6] that the amplitudes for two-body B decays can be 
expressed in terms of 9 diagrams: the color-favored and color-suppressed tree am- 
plitudes T and C , the gluonic-penguin amplitudes P tc and P uc , the color-favored 
and color-suppressed electroweak-penguin (EWP) amplitudes Pew and P%w> the 
annihilation amplitude A, the exchange amplitude E, and the penguin-annihilation 
amplitude PA. These last three all involve the interaction of the spectator quark, 
and are expected to be much smaller than the other diagrams. It is standard to 
neglect them. (Note that the neglect of such diagrams is justified experimentally - 
no annihilation-type or exchange- type decays, such as B^ — > 00, B + — > D s <fi, etc., 
have been observed [TJ.) 

For the three-body decays considered in this paper, we adopt a similar procedure. 
That is, we neglect all annihilation-type diagrams, and express all amplitudes in 
terms of tree, penguin, and EWP diagrams. We assume isospin invariance, but not 
flavor SU(3) symmetry. (It is straightforward to modify our analysis by imposing 
SU(3).) The diagrams are shown in Fig. [TJ A few words of explanation. These 
diagrams are for the decay B — > -r-r-k. There are changes of notation for the other 
decays: 

• For b — > d transitions (B — > KKtt, 7T7T7r), the diagrams are written without 
primes; for b — y s transitions (B — > Knit, KKK), they are written with 
primes. 

• In all diagrams, it is necessary to "pop" a quark pair from the vacuum. It 
is assumed that this pair is uu or dd (= qq); if the popped pair is ss, the 
diagram is written with an additional subscript "s." Thus, for B — > KKn, 
KKK, in the penguin or EWP diagrams with a popped qq pair, the virtual 
particle decays to ss; if the popped quark pair is ss (so the diagram is written 
with an additional subscript "s"), the virtual particle decays to qq. 
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Figure 1: Diagrams contributing to B — > txtxtx. 
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• The subscript "1" indicates that the popped quark pair is between two (non- 
spectator) final-state quarks; the subscript "2" indicates that the popped quark 
pair is between two final-state quarks including the spectator. 

In principle, one can also include the gluonic-penguin diagrams in which the popped 
quark pair is between the pair of quarks produced by the gluon. This corresponds 
to the case where the virtual spin-1 gluon decays to two spin-0 mesons (with relative 
angular momentum 1 = 1). In order to account for the color imbalance, additional 
gluons must be exchanged. Although this can take place at low energy, it will still 
suppress these diagrams somewhat, and so we do not include them here. (Note: 
their inclusion does not change any of our conclusions.) 

One important difference compared to two-body .B-decay diagrams is momentum 
dependence. In two-body decays, in the rest frame of the B, the three-momenta of 
the final-state particles are equal and opposite. One does not have the same type 
of behavior in three-body decays. Although the sum of the three-momenta of the 
final particles is zero, there is no constraint on any individual particle. As such, 
the three-body diagrams are momentum dependent, and this must be taken into 
account whenever the diagrams are used. 

3 Dalitz Plots 

In this section, we review certain aspects of the Dalitz-plot analysis. To illustrate 
these, we focus on the decay B + — > K + 7t~tt + [7J. However, a similar type of analysis 
can be applied to any three-body B decay. 

B + — > K + 7t~tt + can take place via intermediate resonances, as well as non- 
resonant decays. The events in the Dalitz plot are therefore described by the fol- 
lowing two variables: 

X = m K+7r- = (Pk+ +Ptt-) 2 , 

V = m l+n~ = (Pn+ +P7r-f ■ (1) 

Now, one of the great advantages of a Dalitz-plot analysis is that it allows one to 
extract the full amplitude of the decay. To this end, we write 

M(B + K+7T-7T+) = J2c j e i ^F J {x,y) , (2) 

j 

where the sum is over all decay modes (resonant and non- resonant). Cj and 9j are 
the magnitude and phase of the j contribution, respectively, measured relative to 
one of the contributing channels. The distributions Fj, which depend on x and 
y, describe the dynamics of the individual decay amplitudes. In the experimental 
analyses, these take different (known) forms for the various contributions. The key 
point is that a maximum likelihood fit over the entire Dalitz plot gives the best 
values of the Cj and 9j. Thus, the decay amplitude can be obtained. 
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In this paper, the following issue is of central importance. In B + — > K + 7t~tt + , 
since the 7r's are identical particles under isospin, the overall 7r~7r + wavefunction 
must be symmetric. If the tttt pair is in a state of even (odd) isospin, the wavefunc- 
tion (or, equivalently, the B + — > K + h~'k + decay amplitude) must be symmetric 
(antisymmetric) under the exchange p 7T + -H- p n -. Unfortunately, the amplitude of 
Eq. (j2J) does not possess such a symmetry. 

It is the use of the parameters x and y which is problematic. A better choice of 
variables would be s + and s_, where 

s + = m 2 K+7T+ = (p K + +p n +) 2 , 
x = s- = m 2 K+n - = (p K + + p 1T ~f . (3) 

Now, under the exchange p n + <-> p^-, we simply have s + -H- s_. Thus, if we had 
started with the amplitude Ai(B + — > K + h~'k + ) = g(s + , s_), the symmetric combi- 
nation would be 4j[g(s + , s_) +5>(s_, s+)], i.e. it would correspond to the production 

of the 7r~7r + pair with a symmetric wavefunction; ^[g(s + , s_) — gr(s_,s + )] would 
be antisymmetric. 

The problem is that the wavefunction of Eq. (j2J) is not given in terms of s + and 
s_. Fortunately, there is a resolution to this problem: the independent Mandelstam 
variables y, s + and s_ satisfy 

y = m 2 B + 2m^ + - s + - s_ . (4) 

This implies that f(x,y) = /(s_,y) = /(s_, m| + 2m^ +m|-+ -s + -s_) = flr(s+, s_). 
Given the decay amplitude Ai(x, y) of Eq. one can therefore easily construct the 
amplitude which is symmetric/antisymmetric in p n + -h- p n -. The same method ap- 
plies to other B — > K-Kit decays, and indeed to all three-body decays. Thus, if there 
are identical particles in the final state, the 5-decay Dalitz plot allows us to construct 
the amplitude for the production of these particles in a symmetric/antisymmetric 
state. 

Above, we argued that the Dalitz-plot analysis allows one to obtain the amplitude 
M. of any three-body B decay. Actually, this is not quite accurate - the global phase 
of the amplitude is undetermined. Thus, it is really \M\ which should be compared 
with theory. Similarly, one can obtain \M.\ of the CP-conjugate decay. In the 
rest of the paper, we refer to the momentum-dependent branching ratio and direct 
CP asymmetry of a particular decay. These are proportional to |A^| 2 + |A^| 2 and 
lA^I 2 — (A^l 2 , respectively. Finally, for a self-conjugate final state such as K°7t + tt~ 
(where the K° is seen as Kg), the momentum-dependent indirect CP asymmetrjQ 
can be measured, and gives Ai*Ai for this decay. 

4 The indirect CP asymmetry depends on the CP of the final state, and a-priori K°ir + ir~ is a 
mixture of CP + and CP — . However, the separation of symmetric and antisymmetric 7T7t states 
also fixes the final-state CP: K° (iris) sym and K (inr)anti have CP + and — , respectively. 



5 



4 B — > Knn Decays 



We begin with B — > Kirn decays, a b — > s transition. There are six processes: 
B+ -> K+7V+7V-, B + -> IT+ttV, P+ -> T^tt+tt , P^ -> K + n°n-, B a d -> K°n + 7r-, 
B% i^ 7r°7r . In all of these, the overall wavefunction of the final 7T7r pair must 
be symmetrized with respect to the exchange of these two particles. There are two 
possibilities. If the relative angular momentum is even (odd), the isospin state must 
be symmetric (antisymmetric). We refer to these two 

jsym and janti Ag 

shown in Sec. [31 they can be determined experimentally. We discuss them in turn. 

We first consider I^ m , i.e. I = (0,2). The final state has / = §,§, or §. The 
P-meson has I — 2 and the weak Hamiltonian has AI = or 1. The final state 
with I = | cannot be reached. So there are three different ways of getting to the 
final state. Given that there are six decays, this means that there should be three 
relations among their amplitudes. This conclusion is an exact result; the relations 
can be found by applying the Wigner-Eckart theorem: 

A(B + K\ + 7t°) sym = -A(B° d ^K + iT%-) symj (5) 
V2A(B+ K\ + n°) sym = A(B° d K\+ij-) sym + y/2 A(B° d K°7r%°) sym , 
V2A(B° d -> K + u\-) sym = A(B + X + 7r + 7r-) s , m + v^4(B+ -> tf + 7T°7r% m . 

These relations were first given (implicitly) in Ref. [2]. The subscript 'sym' indicates 
that the 7T7t isospin state is symmetrized. 

In terms of diagrams, the amplitudes are given by 



V2A(B+ -+ K°7r + n°) sym = -T[e l1 - C'^ 1 + P' EW2 + P) 
A(B° d ^K°n+n-) sym = - C x e* - P' uc ^ + J* 



IC 

EWl > 



tc 

, 1 p/ I p/C 1 p/C 

~r £ EWL ~ t ~ % EWl ~^ r EW2 > 



I P' -LP' _|_ _ p'C , 1 p/C 
j r £Fl ^ r EW2 * ^ r EWl * <^ r EW2 t 



A{B+ ^ K+K+K-) sym = —T 2 e %1 — Cje 17 — P' uc e 11 + P[ c 

-P' --P ,c + - 

^ r EWl ^ r EWl * g 

V2A(B + -»• K + 7i°7T°) sym = T[e^ + + C[e^ + C 2 e^ + P' uc e^ - P' t( 



I P' P'C _l_ P /C< 

r ~- r EWl o EWl o EW2 > 



fc 



1 p/ p/ 2 p/C 2 p/C 

g- r BW r l - r EVK2 ^ EWl ~^ r EW2 > 



>/2A(SS ir + 7r 7r-) s , m = r i ^ + C' 2 ^-P' EW2 -P'£ wl , (6) 

where P' = P[ + P^. (Note: all amplitudes have been multiplied by a/2.) Above we 
have explicitly written the weak-phase dependence (including the minus sign from 



6 



V t * b Vt s [P/ c and EWP's]), while the diagrams contain strong phases. (The phase infor- 
mation in the Cabibbo-Kobayashi-Maskawa quark mixing matrix is conventionally 
parametrized in terms of the unitarity triangle, in which the interior (CP- violating) 
angles are known as a, (3 and 7 [8 J.) It is straightforward to verify that the three 
relations of Eq. (jSJ) are reproduced. Thus, in this case, there is no difference between 
the exact and diagrammatic amplitude relations. 

We now turn to I%% tl , i.e. 1 = 1. Here there are four processes: B + — > K + 7i + 7i^ , 
B + — > K 7t + 7t°, B d — > K + 7t°7t~ , B d — > K°7t + 7t~ (one cannot antisymmetrize a 7r°7r° 
state). The final state has I = | or |, so there are still three different paths to 
get to the final state. We therefore expect one relation among the four amplitudes. 
Ref. [2] notes that it is similar to that in B — > nK: 

V2A(B + -> K + n + n-) antl + A(B + -> K\ + n°) anti = 

V2A(B° d -+ K\ + n-) anti + A(B° d -> K + Tt\-) antl , (7) 

where the subscript l antV indicates that the tttt isospin state is antisymmetrized. 

Writing the amplitudes in terms of diagrams is a bit more complicated because 
antisymmetrization is involved. Depending on the order of the pions, there might 
be an extra minus sign. To account for this, we use the following prescription: 

• All diagrams with the pions in order of decreasing charge from top to bottom 
are unmodified; all diagrams with the pions in order of increasing charge from 
top to bottom get an additional factor of —1. 

This requires that diagrams always be drawn the same way. For example, the 
spectator quark for all tree diagrams should always appear in the same place (e.g. 
at the bottom of the diagram), and the decay products of the neutral bosons in 
penguin and EWP diagrams should always appear in the same order (e.g. quark on 
top, antiquark on the bottom). 

With this rule, the amplitudes take the formEl 



V2A(B + -> K\ + 7r°) anti = -T[^ - C' 2 e^ - 2P^ + 2P U 

1 c 2 



tc 



Pi/ p'C 1 p'C 

EW2 ~ o^EWl t o' r EW2 > 



A(B° d ^ K\ + n-) antt = -T[^ -Cie* - P^e* + P u 



tc 

, p> 2 p'C , 1 p'C 

1* r EWl ^ r EWl ^ jr EW2 ■> 



A(B + -+ K + n + n-) anti = T 3 e* - C[e^ + P' uc e^ - P' ti 



tc 



1 pi 1 p'C 1 2 p'C 

' r EWl ~^ r EWl ' ~^ r EW2 1 



5 Note: even though the diagrams of Eq. ([5]) have the same names as those of Eq. ©, they are 
not the same diagrams. That is, in general, they take different values. 
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V2A(B° d -> K + 7i\-) anti = T[^< + 2T^ - C 2 e* + 2P' uc e l ~< - 2P' tc 

~ P'e\V2 + qP'ewi + ~^Pew2 ■ (8) 

(As above, all amplitudes have been multiplied by a/2.) The relation of Eq. ([7]) is 
reproduced. Therefore, there is no difference between the exact and diagrammatic 
amplitude relations in the antisymmetric case. 

4.1 Resonances 

It is possible that the B decays to an intermediate on-shell M\M 2 state, which then 
subsequently decays to Kirn. Examples of such resonances are M\M 2 = Kp, K*ir, 
Kf (980). The question now is: how does the diagrammatic analysis presented 
above jibe with resonant decays? To answer this, we examine the resonances in 
turn. 

Consider first M X M 2 = Kp. The four decays are B + K + p°, B + K°p + , 
B d — >• K°p°, B d — > K + p~ , whose amplitudes take the form 

V2A(B + K+p°) = -Tyf* - P <* - P' ucy ^ + P[ cy + P' EmP + 2 -P'g wy , 

A(B + ^K°p + ) = P' ucy e^ - P' tcy + \p'§ wy , 
s/2A(B° d -> K°p°) = -C' P ^ + P^yt* - P[ cy + P' EWjP + ~P' E c wy , 

A(B° d ^K + p-) = -TW~< - Pi cy e^ + P> tcy + 2 -P>g wy , (9) 

where the subscript P or V indicates which final-state meson [pseudoscalar (K) or 
vector (p)} contains the spectator quark of the B meson [9]. (Note that the diagrams 
which describe resonant decays are a subset of those used for B — > Kixix (Fig. [Q. 
Above, the diagram Dy (Dp) is the same as D 2 (-Di)-) The relation among the 
amplitudes is 

V2A(B + -> K + p°) + A(B + K°p + ) = 

V2A(B° d -> K°p°) + A(B d -> K + p-) . (10) 

Given that p° — > 7r + 7r~, p + — > 7t + tt° and p~ — > 7r°7r~, this reproduces Eq. ([7]), which 
is the relation for the antisymmetric tttt isospin state. This makes sense, since the 
p decays to {mT) anti . 

Consider now M X M 2 = Kf (980). There are two decays: B + K + f (980) and 
B d — > K° / (980) . It is straightforward to show that there is no relation between the 
two amplitudes. However, the /o(980) decays to a pion pair in a symmetric isospin 
state, with A(f — > (n + n^) sym ) = —\/2A(f — > 7r°7r°). This leads to 

A{B° d K\ + 7i-) + V2A(B° d -> K\°n°) = , 
A(B + ir + 7r + 7r-) + ^^(5+ K + ir\°) = . (11) 
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Given that the if /o (980) resonance does not contribute to A(B + — > K°tt + tt°), 
A(B° d ->■ K + n°7i-) or A(P+ -»■ if°7r+7r ), the decays B -»■ if/o(980) -»■ ifvrvr 
satisfy Eq. flSJ), which are the relations for the symmetric tttt isospin state. 

Finally, consider M 1 M 2 = K*tx. The four decays are B + — > K*°tt + , B + — > 
K*+tt°, B° d K* + 7i~, B° d K*°tt°. The amplitudes are [9] 

A{B + ^K*\ + ) = PL, P e l " - K,P + Ip'e C w,p , 
V2A{B+ if *+tt°) = -r p e« - C' v ^ - P' ucP e^ + P[ cP + P' Ewy + %P'^ P , 



3 



A(B° d -> iT+vr-) = -T^ - P^ P e^ + P; jF + ^P^ )P , 



V^A(BS iTV) = -C(,e*r + P^ - P' tCjP + P' EWy + \p'^ P . (12) 

The relation among the amplitudes is 

A(B + ->■ K*°n + ) + V2A(B + ->• K* + n°) = 

A(B d ->■ ir+vr-) + V2A(P° -). iTV) . (13) 

Now, the if* decays to if7r, and both charge assignments are allowed: 



if*+ -> Jl/3if + 7T U - J2/3 if °7T 



1 



if* u ->• y/2/3 K+ic- - ^1/3K%° . (14) 

There are therefore several if *7r contributions to a particular Ktttt final state. How- 
ever, one never reproduces the relations in Eqs. (|5]) or (JTJ). This reflects the fact that 
this resonance contributes to both (7nr) sym and (mr) ant i. 

Still, it is instructive to examine the relation obtained when the resonance decays. 
This is obtained by inserting Eq. (I14p into Eq. (fl3|) . When the tttt pair is in a 
symmetric isospin state, one has 

V2A(B + -»■ if+TT+vr-) - 3A{B + K°tt + tt°) + V2A(B + K + tt°tt°) = 
3A(B° d -> if + 7T°7r-) - v/2A(P° -> if°7T + 7r-) - ^(5° -> if°7r°7r ) .(15) 

This is obviously not the same as Eq. (j5]). This is because there are only four 
£> — > K*tt decays (and not six, as in B — > Ktttt), and so there is only one relation 
among the if tttt decays. 

On the other hand, the case where the tttt pair is in an antisymmetric isospin state 
is more interesting. For I%™ u , amplitudes to final states with two 7r°'s are zero. Also, 
there is an additional factor of —1 if the pions are in order of increasing charge from 
top to bottom. Taking the if* in B — >■ K*tt to be on top of the tt, the amplitudes 
A(B + -> K + ir + ir-) K . 0v+ , A(B + -»■ if°7r + 7r°)^o 7r+ and A(B° d -»■ K^tt°tt~) k ^ all 
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get an extra minus sign (the subscript indicates the resonance which gives rise to 
the final state). When these are taken into account, the insertion of Eq. (j!4p into 
Eq. (Tl3l) gives the relation in Eq. ©. We therefore see that the B — > Knii amplitude 
relation is reproduced by B — > K*n decays for the /"" tl case. 

The point here is that it is useful to consider the entire B — > M\M.2 — > Ktttt 
decay chain, and that the distinction between and is important, even for 
resonances. 

4.2 Penguin Dominance 

In general, the dominant contribution to b — > s transitions comes from the penguin 
amplitude. In Ref. [I], Gronau and Rosner explore the consequences for B — > Kmr 
decays of assuming penguin dominance and neglecting all other contributions. They 
note that, in this limit, the amplitudes must respect isospin reflection (i.e. u -H- d), 
which implies that 

A{B + -> K + 7T + 7T-) = A(B° d K°7T + 7l-) , 

A(B + ->> K\ + tt°) = A(B° d K + tt\-) , 
A(B° d -> K 7r°n°) = A(B + -> /^ + 7r 7r ) , (16) 

up to possible relative signs. They find that, on the whole, the data respect these 
relations. 

The expression of the amplitudes in terms of diagrams allows us to go beyond 
these results. Using the method of Sec. [3] to distinguish I^ 71 and I%™ u , it is possible 
to consider the two cases separately, under the condition that only the diagram P( c 
is retained in the amplitudes. 

In the symmetric scenario, we have the following predictions: 

A(B + -> K\ + tt°) = A(B° d -> K + ti\-) = , 
A(B + ->■ K+7T+7T-) = A(B d -> K°7i + 7i-) 

= - V2A(B° d -> K°7T 7r ) = -V2A(B + K + ir°n°) . (17) 

And in the antisymmetric scenario, we have 

A(B° d ^ K°n\°) = A(B+^K+n\°) = 0, 
A(B + -)• K\ + tt°) = —A(B d -> K + tt\-) 

= - V2A{B + -> K + ti + ti-) = V2A(B° d -> ^tt+tt") . (18) 

These provide further tests of the SM. 

In fact, several of these decays have been measured: B + — > K + n + 7i~ [7J [TO] . 
B° d -»■ i^vr+vr- [HJ, and 5g -> K+tt^-k- [12J. We can therefore test some of the 
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above relations. Specifically, in terms of branching ratios (integrated over the entire 
Dalitz plot), the predictions are 



B{K+Ti\-) sym = 0, 

B(K+n + n-) sym = (r + /r )B(K\ + n-) sym , 
^(T + /T )B(K + n\-) anti = B(K + n + n-) anU = (T + /T )B(K\ + n-) anU . (19) 

We determine the symmetric and antisymmetric amplitudes for the three decays 
using the Dalitz-plot method described in Sec. El Consider first B + — > K + 7t + tt~ . 
We write this amplitude in terms of x = (pk+ + P-K+) 2 an d y = (pk+ + P-k-) 2 - 
Given the decay amplitude f(x,y), the symmetric amplitude is taken to be f sym = 
■^(f{x, y) + f(y, x)), and we compute the integral of \f sy m\ 2 and \f\ 2 over the Dalitz 
plolE|. A similar procedure is carried out for the antisymmetric amplitude f an u = 
-?=(/(:r, y) — f(y,x)). The other two decays are treated in the same way. 

Although the full amplitudes for B + — > K + -k + 'k~ and B% — > iC 7r + 7r _ are split 
roughly equally between symmetric and antisymmetric, the same is not true for 

T(K + n + n-) sym = O.Q5T(K + n + n-) , 
T(K°7r + 7i-) sym = 0.68r(K°7r+7r-) , 

r(^ + 7r°7r-) s , m = 0.11 T(K + n\-) . (20) 

With these, we obtain 

B(K + n°n-) sym = (4.0 ±0.3) x 10~ 6 , 
B(K + n + n-) sym = (33.3 ± 2.0) x 10 



-6 



r + /T )B{K\ + n-) sym = (36.4 ± 1.5) x KT 6 

( 

2 



I (t + /t ) B(K + n\-) antl = (17.1 ± 1.3) x 10 



B(K + 7r + n-) antl = (17.6 ± 1.0) x 10- 6 
{r + /r )B(K°n + n-) antl = (17.0 ± 0.7) x KT 6 



(21) 



(Note that the above errors do not include the errors in the parameters obtained from 
the Dalitz-plot analyses of the three decays.) We therefore see that the data agree 
with the predictions of Eq. ( 1T9"|) . In particular, B(K + 7i n^) sym is indeed greatly 
suppressed, in agreement with the SM. 



6 Note that, because of the coefficient ^ in f sym , one must integrate over only half of the 

Dalitz plot to avoid double counting. Alternatively, f sym can be defined with a factor i, and one 
integrates over the entire Dalitz plot. There are no such issues with /. 
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4.3 Weak-Phase Information 



Since the expressions for the decay amplitudes include the weak phase 7, it is natural 
to ask whether 7 can be extracted from measurements of B — >■ Ktctc decays. The 
answer is 'yes' if the number of unknown theoretical parameters in the amplitudes 
is less than or equal to the number of observables. In performing this comparison, 
we examine separately the and scenarios. 

Consider first the 1^™ case. Here there are six B — > Kmr decays. On the other 
hand, the first relation in Eq. flS]) shows that the amplitudes for B + — > K 7t + 7t 
and B d — > K + tt ti~ are equal (up to a sign), so that there are only five independent 
decays. The Dalitz-plot analyses of these decays allow one to obtain the momentum- 
dependent branching ratios and direct CP asymmetries of B + — > K + tt + tt~, B + — > 
JT+ttV, B Q d -> K + 7T°7r-, B a d ->• K°7i + n-, and B° d K°n 7i . In addition, one can 
measure the momentum-dependent indirect CP asymmetry of B d — > K°7t + tt~ . (The 
indirect CP asymmetry of B d — > K°tt°it will be very difficult, if not impossible, to 
measure.) Thus, there are essentially 11 (momentum-dependent) observables in I^" 1 
B — > Kim decays. 

For the case of there are four decays, yielding 9 observables: the momentum- 
dependent branching ratios and direct CP asymmetries of B + — > K + 7t + tt~ , B + — > 
fr°7T + 7r , B d — > K + 7T°7T~ , B d — > K°7t + 7t~, and the momentum-dependent indirect 
CP asymmetry of B d — > K 7t + tt~ . Since this is fewer than above, we conclude that 
the I^ 71 scenario is the more promising for extracting 7. 

The six i^f™ amplitudes are given in Eq. fl6]). Although there are a large number 
of diagrams in these amplitudes, they can be combined into a smaller number of 
effective diagrams: 



V2A{B + ->• K°7i + 7i Q ) 
A(B° d ->• K\ + n-) 
V2A(B° d -> K\\°) 



sym 



sym 



sym 



—T'.P^ 4- P'pi'y — P' -L P' _|_ p' 

_ p* Jn _i_ P' _ P' 

r aF + r b r EW,a 1 

t'p^ 4- t'p^ 4- p'p ir y — p' — p' 

T'p^ 4- T'p in l — P' — P' 
1 a t + 1 b e r EW,a r EW,b j 



_rpi vy _ rpl 17 1 pi 1 pi 

1 a fi 1 6 e ^ r EW,a ^ r R 



A(5 + -)■ K+7T + 7T") 
>/2A(B + ->• ^ + 7T 7T ) 
V2A(B° d -)• ir + 7T 7T-) 



sym 



sym 



sym 



(22) 



where 



P' 




(23) 
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The amplitudes can therefore be written in terms of 6 effective diagrams. This 
corresponds to 12 theoretical parameter^]: 6 magnitudes of diagrams, 5 relative 
(strong) phases, and 7. We remind the reader that the diagrams are momentum 
dependent. This does not pose a problem. They will be determined via a fit to the 
data. But since the experimental observables are themselves momentum dependent, 
the fit will yield the momentum dependence of each diagram. 

Unfortunately, as noted above, there are only 11 experimental observables. There- 
fore, in order to extract weak-phase information (7), one requires additional input. 

A previous analysis made an attempt in this direction. In 2003, Deshpande, 
Sinha and Sinha (DSS) wrote schematic expressions for the symmetric B — > Kw 
amplitudes, including tree and EWP contributions [13] . Now, in B — > ttK decays, it 
was shown that, under flavor SU(3) symmetry, the EWP diagrams are proportional 
to the tree diagrams (apart from their weak phases) [H]. DSS assumed that the 
EWP and tree contributions to B + — > fC°7r + 7r° are related in the same way. This 
gives the additional input, and allows the measurement of 7. Unfortunately, it was 
subsequently noted that the assumed EWP-tree relation in Knn does not hold [TJ] , 
so that 7 cannot be extracted. This is the present situation. 

In fact, the situation can be remedied. Referring to the B® — > K 7i + 7T° amplitude 
in Eq. DSS made the assumption that T{ + C' 2 is related to P' EW2 + P'ewii an< ^ 
this was shown not to be true. We agree with this. However, there are other EWP- 
tree relations which do hold, and their inclusion does allow the extraction of 7. The 
full derivation is rather complicated, and so we present this in a separate paper [16J. 

Finally, we note that there is another method for obtaining 7 from B — > Kmr 
decays. In two-body b — )■ s B decays, the diagrams are expected to obey the 
approximate hierarchy [B] 

1 P' 

A : T , P EW , 

A 2 : Pgr, (24) 

where A ~ 0.2. If the three-body decay diagrams obey a similar hierarchy, one can 
neglect C[, C 2 , P' uc1 P'ewii P'ew2i an d incur only a ~ 5% theoretical error. But if 
these diagrams are neglected, then two of the effective diagrams vanish: P' EWa —> 
and T' h — P' a — > [Eq. (j2"3"|) ]. In this case, the amplitudes can be written in terms 
of 4 effective diagrams, corresponding to 8 theoretical parameters: 4 magnitudes of 
diagrams, 3 relative (strong) phases, and 7. Given that there are 11 experimental 
observables, the weak phase 7 can be extractecH. 

7 In fact, there is another theoretical parameter - the phase of Bj-Bj mixing, /3, enters in the 
expression for the indirect CP asymmetry. However, the value for j3 can be taken from the indirect 
CP asymmetry in B° d -> J/^K S [5]- 

8 This technique does not work when the nir pair is in an antisymmetric state of isospin. In this 
case, there are still more theoretical unknowns than observables, so that 7 cannot be extracted. 
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The downside of this method is that it is difficult to test the assumption that 
certain diagrams are negligible. Indeed, the presence of resonances may change 
the hierarchy. In light of this, the theoretical error is uncertain, and this must be 
addressed if this method is used. 



5 B KKK Decays 

We now turn to B — > KKK decays, also a b — > s transition. The four processes 
are: B + -> K + K + K~, B+ -> K + K°K°, B° d -> K + K°K~, B° d -> K°K°K°. Here 
the overall wavefunction of the final KK pair must be symmetrized. If the relative 
angular momentum is even, the isospin state must be symmetric {I — 1); if it is 
odd, the isospin state must be antisymmetric (1 = 0). 

For the symmetric case, the final state has I = \ or |, so there are three different 
ways of reaching it. There should therefore be one relation among the four decay 
amplitudes. From the Wigner-Eckart theorem, it is 

A(B+ -> K + K + K-) sym + V2A(B + -> K + K°K°) sym = 

V2A(B° d -> K + K°K-) sym + A(B° d K°K K°) sym . (25) 

In terms of diagrams, the amplitudes are given by 

A(B + ^K + K + K-) sym = -T>^ - - P' uc ji + P' tc 

, 2 P / 1 P / , 2 p/C 1 p/C 

"T g-'BVKl.S g-'EH/l jj r W2,S g^EWl 5 

>/2A(B + -> K + K°K°) sym = P' uc ^-P[ c 

1111 
i _ P' J._P' _l _ pc , i p/C 

>/2A(5S K+K°K-) sym = -T^-C^-P^ + P/, (26) 



, 2 p/ 1 p/ , 2 p/C 1 p/C 

"r ^ r EWl,s c^ r EWl ^ r EW2,s r^ r EWl > 



A(B° K K ^°) SJ/m = P' uc e-'-P[ c 

1111 

_L _ P/ J._P' I _ P'C _L _ P'C 

where P' = P' 2 s + P{ . It is straightforward to verify that the relation of Eq. ( I25l) is 
reproduced. On the other hand, one sees that there are, in fact, two relations: 

A(B + -> K + K + K-) sym = V2A(B° d ->• K + K°K~) sym , 

V2A(B + ^K + K°K°) sym = A(B d — >■ K K K ) sym . (27) 

What's happening is the following. Eq. f )25|) is exact. However, when annihilation- 
type diagrams are neglected - as is done in our diagrammatic expressions of ampli- 
tudes - then one finds the two relations above. This is an example of how one can 
go beyond the exact relations if certain negligible diagrams are dropped. 



14 



In order to test these relations, it is necessary to isolate the symmetric piece 
of the decay amplitudes. B + — > K + K + K~ and B d — > K°K°K° are automatically 
symmetric since the final states contain truly identical particles. On the other hand, 
for B d — > K + K°K~ and B + — > K + K°K°, the symmetric amplitude can be obtained 
using the Dalitz-plot method of Sec. |3j Now, the Dalitz plot of B d — > K + K°K~ has 
already been measured [T71 [18]. This allows us to test the first relation in Eq. ff27|) . 

We use the Dalitz-plot analysis of B d — > K + K$K~ given in Ref. [17], with 
A(B° d K + K°K~) = V2A(B° d K + K S K~). We find T(B° d K+K°K-) sym = 
0.57T(B° d K + K°K~). This then gives 

2 (r+/r ) B(B° d ->• K + K°K-) sym = (30.0 ± 2.8) x 10" 6 . (28) 

(Note that the above error does not include the errors in the parameters obtained 
from the Dalitz-plot analysis of Ref. [UJO This is to be compared with [1] 

B(B + -> K + K + K~) = (32.5 ± 1.5) x 10~ 6 . (29) 

We therefore see that the first relation in Eq. (|27j) is satisfied. This supports our 
assumption that annihilation-type diagrams are negligible. 

In the antisymmetric case, there are only two decays: B + — > K + K°K° and 
B° d K + K°K~. A(B + -> K+K+K-) and A(B° d -» K°K°K°) vanish because 
there is no way of antisymmetrizing the K + K + or K°K° pair. Here the final state 
has I = |, and there are two different ways of reaching it. We therefore expect no 
relation between the amplitudes. 

In order to write the amplitudes in terms of diagrams, we have to antisymmetrize 
the K + -K° state. As was done for Kirir, we adopt the following rule: all diagrams 
with the K + -K° in order of decreasing charge from top to bottom are unmodified; 
all diagrams with the K + -K° in order of increasing charge from top to bottom get 
an additional factor of —1. The amplitudes (multiplied by a/2) are then given by 

V2A(B+ K + K°K°) anU = -P' uc e l -< + P' tc 

1 P' _ 1 P' _i_ 1 P'C I 1 p'C 

c^ r EWl,s <^ r EWl <^ r EW2,s < r^ r EWl ' 

V2A(B° d ^K + K°K-) anti = -7*,e*r + C^e* - P' uc e^ + P' tc (30) 

, 2 p> _ 1 p' _ 2 P'C i 1 P'C 

i g-'EVKl.s ~^ r EW\ ~^ r EW1,s ' ~^ r EWl ' 

As expected, there is no relation between these two amplitudes. 
5.1 Penguin Dominance 

Assuming penguin dominance, Gronau and Rosner find that isospin reflection im- 
plies the following equalities jl]: 

A(B + -)• K + K + K-) = —A(B d -)• K°K°K°) , 
A(5+ ->■ K + K°K°) = —A(B d K + K°K-) . (31) 
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By distinguishing the symmetric and antisymmetric isospin states, it is possible to 
go beyond these predictions. In the symmetric scenario, if only P' tc is retained, we 
predict 

A(B + -> K+K+K-) = -A(B° d -> K°K°K°) 

= - V2A(B + -»■ K + K°K°) = V2A(B° d K + K°K~) . (32) 

(Note: the relations given in Eq. ( 127]) actually hold for all diagrams, not just P( c .) 
As discussed above, the present data confirm the relation A(B + — > K + K + K~) = 
\/2A(B d — > K + K°K~). In the antisymmetric scenario, we have only A(B + — > 
K + K°K°) = A(B° d ->■ K + K°K~). As with Ktctt decays, these provide further tests 
of the SM which. 



5.2 Isospin Amplitudes 

In Ref. [3], Gronau and Rosner (GR) write the amplitudes for B — > KKK decays in 
terms of isospin amplitudes. It is instructive to compare this with the diagrammatic 
description. 

As described above, there are five independent isospin amplitudes, denoted 

by A^f K),If = (I(KK),I f \Al\±), where I(KK) is the isospin of the KK pair 

[I(KK) = 1 (0) is symmetric (antisymmetric)], // is the isospin of the final state, 

o,i ii o,i 

and the weak Hamiltonian has AI = or 1. They are listed as A ' 2 , A ' 2 , A{ 2 , 

1 I 1 3 

a\ 2 ,a\ 2 . 

As noted by GR, the B — > KKK amplitudes depend on the kaons' momenta. 
The amplitudes for B + — > K + K°K° and B d — > K + K°K~ take different values when 
the K + and K° momenta are exchanged. Thus, GR obtain expressions for six decay 
amplitudes in terms of the five isospin amplitudes: 

A(B + ^K + K + K-) PlP2m = 2AY 2 - 2AY 2 + A 1 / 2 , 

A(B° d ^ K K K°) piP2P3 = -2A 1 / 2 - 2A^ 2 + A 1 / 2 , 

A(B + -> K + K°K°) PlP2P3 = A / 2 - AY 2 - A Q ^ + a\'~ 2 + A^ , 

A{B + -> K + K°K°) P2PlP3 = -A / 2 - A 1 / 2 + A°/* + A 1 / 2 + A 1 / 2 , 

A(B° d -> K + K°K-) P1P2P3 = A / 2 + A 1 / 2 + A?* + A 1 / 2 + A 1 / 2 , 

A(B° d K + K°K-) P2P1P3 = -A / 2 + A 1 / 2 - A*'" 2 + A^ + A 1 ,' 1 . (33) 

The above amplitudes are related to those of Eqs. ( 1261) and ( 1301) as follows: 

A{B + ^K+K + K-) sym = A(B+ ^K + K + K-) P1P2P3 , 
A(B d ^K K K°) sym = A(B° d ^K K K°) piP2P3 , 
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V2A(B + -> K + K°K^ 



I sym 



A(B + -> K + K K Q ) piP2P3 + A(B + -»■ K + K°K°] 



P2P1P3 ' 



V2A(B° d ^K + K°K~) sym = 

A(B° d -> K + K°K~ 



IP\P2PZ 



/P2P1P3 ' 



\/2A(fi + -> K + K°K^ 



lanti 



A(B + K + K°K°) PlP2P3 - A{B + -4 K + K°K°) P2Pip , A , 
V2A(B d ^K + K°K~) antl = 

A(B° d T^ ^)™ - -> ^ + K°^-) P2P1P3 . (34) 

Now, because there are six decay amplitudes, but only five isospin amplitudes, 
there must be a relation between the decay amplitudes. GR give this relation as 

-»■ K + K + K~ 



)piP2P3 A(B v K K K ) PlP2P3 
+ -> K+K°K°] 

A(B° d -> K°K°K\ 1P2P3 + A(B d K+K°K-) PlP2P3 

+ A(5^ -> K + K°K~ 



P2P1P3 



!P2P\P-A 



3A^ . (35) 



This is the same as the relation in Eq. ( )25|) . However, when one expresses the 
amplitudes in terms of diagrams, there are, in fact, two relations instead of one 
[Eq. (jTTJ) ]. This implies that 

i I 1 l 2 

^l' 2 = -^l' 2 , (36) 

so that there are really four independent isospin amplitudes instead of five. As 
described above, the extra relation is a consequence of neglecting the annihilation- 
type diagrams. In other words, the above relation among isospin amplitudes is 
a good approximation, and could not have been deduced without performing a 
diagrammatic analysis. 

It is straightforward to express the remaining isospin amplitudes in terms of 
diagrams: 



^0 



1 2 



An 



a** 



- V - <V - 2p Le 11 + 2 ^ c 

1 9 1 

4- - P' --P' +-P' C 

1 2,s c u l,s e ^ EWX,s ' EW2 



2 p/C 



£VK1 



-T' e *r + C* e* - 2P^ + 2P' tt 



4. 1 p> --P' 

' ^ r EWl,s ^ E W1 



tc 

1 P'C _L 2 p'C 

g- r EW2,s t --Lewi 



-T' p* 7 -i- r" p* 7 4- P' — P' c 
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(37) 



(Recall that, despite their having the same name, the diagrams which contribute to 

1/1.2} o± 

the ^ 2 and A^ 2 ^ isospin amplitudes are not the same - they can have different 

i,- l, 1 

sizes.) In the limit of penguin dominance, A{ 2 and A ' 2 vanish. This is consistent 
with what is found in the previous subsection. 

5.3 Weak-Phase Information 

As was the case for B — y K'ktt decays, the amplitudes contain the weak phase 7, 
and so one wonders if it can be measured in B — y KKK decays. Here the answer 
is 'perhaps'. 

When the isospin state of the KK pair is symmetric, there are four decays. 
However, due to the equality relations in Eq. (I27|) . two of these have the same 
amplitudes as the other two. There are therefore 6 observables: the momentum- 
dependent branching ratios, direct CP asymmetries and indirect CP asymmetries 
of of B® — y K + K°K~ and B® — y K°K°K°. In the antisymmetric scenario, there 
are 5 observables: the momentum-dependent branching ratios and direct CP asym- 
metries of B + — y K + K°K° and B® — y K + K°K~ , and the momentum-dependent 
indirect CP asymmetry of B% — y K + K°K~ . (As with B —y Kirn, the separation 
of symmetric and antisymmetric KK states fixes the CP of the final state for the 
indirect CP asymmetries.) 

However, in either case, the amplitudes [Eqs. (j26|) and ( 130]) ] are written in terms 
of 4 effective diagrams, corresponding to 8 theoretical parameters: 4 magnitudes 
of diagrams, 3 relative (strong) phases, and 7. This is larger than the number of 
observables, and so the weak phase 7 cannot be extracted from B — y KKK decays. 

The best that one can do is to assume the hierarchy of Eq. (I24p . and neglect 
all C, P' uc and P'ew diagrams. This reduces the number of effective diagrams to 
3, which corresponds to 6 theoretical parameters. This is equal to the number of 
observables in the symmetric case, so that 7 can be extracted here, albeit with 
discrete ambiguities. And, as described above, the theoretical error is uncertain. 

6 B — » KKn Decays 

We now consider B —y KKn decays, which are b — y d transitions. Here there 
are seven processes: B + -> K + K--k\ B + ->■ K + K°n°, B + -y K K n + , B° d ->■ 
K + K-n°, 5° -> K + K°n-, 5° -> K°K°7r°, B° -> K°K-n + . There are no identical 
particles in the final state, so here we do not have to distinguish symmetric and 
antisymmetric isospin states. 

In B — y KKtt, the final state has I = 0, I = 1 (twice) or I — 2. The weak 
Hamiltonian has AI = | or |, so there are six paths to the final state. This implies 
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that there is one relation among the seven decay amplitudes. It is 

V2A(B° d -)■ K+K'tt ) + A(B° d ->■ K°K-n + ) - A(B + -> K + K-n + ) 
+ V2A(B° ->■ K°K°tt°) + A(B° d ->■ K+K\~) 
- A(P+ -> K°K\ + ) - V2A(B + -> K + K°ir°) = . (3* 

In terms of diagrams, the amplitudes are given by 

-+ K+K-TI+) = [T 2tS + C hs + P a]UC }e- ia 

1 2 1 c 2 c 



= [T 1)S + C 2iS -P o;uc + P 6;uc ]e- 



c c c 

+ -Pa;tc — P&;tc — P,EW2,s — g-^EWl ~ g °.EWl,s + g °£W2 — g°£W2,s ) 
A(P+ -> K°i? 7T + ) = 

+ -Pb;tc — -^Pewi — -^Pewi,s — -^Pewi,s ~ ^Pew2 j 
V2A(B° d -)■ K + K-n°) = C^e'™ + ^P EW1 -^P EW1 , S , 

A(B° d ^K + K°n-) = [T ljS + P b , uc ] e~* a - P b , tc - ^ W1 , s + \p C EW 2 , 
V2A(B° d ^K K ir°) = [C 2 , fl - P a]UC - P b;uc ] e~ ta (39) 

+ Pa;tc + Pb;tc ~ ~PeW1 ~ ^PeW1,s ~ PeW2,s 

1 P C 1 pC 1 pC 1 pC 

— g-'EU'l g-TEWl.s ~^ r EW2 ^ EW2,S J 

A(P° if°#-7r+) = [T 2)S + P a;uc ] - P a;tc + ^P£ m - j^JW , 

where P a = Pi + P 2iS , p, = Pi |S + P 2 , and all amplitudes have been multiplied by 
e*^. With these expressions, the relation of Eq. ( 1381) is reproduced. 
However, there are, in fact, two relations: 

V2A(B° d -> K + K-7T°) + A(B d -> K°K-tt + ) = A(B + -> K + R-n + ) , 
V^A(P° -)• K°K°n ) + A(B d K + K\-) 

= A(B + ->■ K°K\ + ) + v^A(S + ->• K + K\°) . (40) 

As was the case in P — >■ KKK decays, the (justified) neglect of certain annihilation- 
type diagrams breaks the relation in Eq. ( 155]) into two. 
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6.1 T Dominance 



In two-body B decays, T is the dominant diagram in b — > d transitions. Assuming 
this also holds in three-body B decays, we have the following predictions: 

A(B + -> K + K-n + ) = A{B° d -> K°K-n + ) , 
V2A(B + K + K°tt°) = A(B° d ->■ K + K°n~) , 
-> K°K°tt + ) = A(B° d -> K + K-ir°) = A(B° d -> K°K\°) ~ . (41) 

These are tests of the SM which can be carried out once these decays are measured. 



6.2 Weak- Phase Information 



There are seven S — >■ KKir decays, which yield 16 observables: the branching ratios 
and direct CP asymmetries of 5+ -> K + K~n + , B + -> K + K°n , B + -> K K°n + , 
B° d ->■ K+K-n°, 5° -> K + K°tt-, B° d -> K K°tt°, B° d -> K°K-ir+, and the indirect 
CP asymmetries of 5° -> K + K~n°, 5° ->■ K°K n . 

The S —7- KK-k amplitudes in Eq. ( 13 9 p can be written in terms of 10 effective 
diagrams: 



V2A(B + -> K+ifV 

A(5° ->■ K + K\- 
V2A(B° d i^ifV 
A(5° -> K°K-n + 



-7T 



[-D 5 + £> r + D 9 ] e - fa -£> 6 



10 



(42) 



where 





= Cl.s ) 






= g-TBWl _ g"Wl,S 




£ 3 


= -^2,s "I - Pa;uc j 




£4 


- P _L 1 pC 

— -t a; tc ~T~ g-T^vyi 


V 

2 -'£11/2, s ) 


£5 


= -^l,s "I - Pb;uc 1 






— P _L 1 pC 

— -r&jic -r g jr EW2 


- P c 




= Pb;uc 1 






P l P 1 


P 1 

"Wl,s — g 
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£W/2 



2 r BWl,« > 



Dg = Ti jS + C^s — P a ;uc + Pb;uc , 
DlO = -Pajte — Pb;tc — PEW2,s ~ g^Wl ~ ^PeW1,s + ^EW2 ~ ^PeW2,s ■ (43) 

This corresponds to 20 theoretical parameters: 10 magnitudes of diagrams, 9 relative 
(strong) phases, and a. With only 16 observables, a cannot be extracted. 

We therefore need additional input. Fortunately, we have some, similar to that in 
Sees. 14.31 and 15.31 In two-body b — >■ d B decays, the diagrams obey the approximate 
hierarchy [6] 



1 




A 


n p p 

^ ; 1 tC) 1 UC 


A 2 


Pew , 


A 3 


pC 

r EW ■ 



If the three-body decay diagrams obey a similar hierarchy, all EWP diagrams can 
be neglected, leading to an error of only ~ 5%. In this limit, we have D2 = 0, 
D 8 = —Dq, and D w = — _D 4 + D e . So the number of independent diagrams is 
reduced to 7, i.e. 14 theoretical parameters^. Thus, by measuring the observables 
in B — > KKtt decays, weak-phase information can be obtained. In fact, not all 
16 observables are necessary. Experimentally, this is not easy, but it is at least 
theoretically possible. Of course, as in Sees. 14.31 and 15. 3[ the theoretical error is 
uncertain, since it is difficult to test the hierarchy of diagrams. 



7 B nnn Decays 

Finally, we examine B — > it-kit decays, also a b — > d transition. There are four 
processes: B® -> 7r 7r°7r°, B + -> 7r + 7r°7r°, B + -> it~it + it + , B Q d -> it + it q it~ . In 
contrast to the other decays, here the final state includes three identical particles 
under isospin, so that the six permutations of these particles (the group S3) must 
be considered. Numbering the particles 1, 2, 3, the six possible orders are 123, 132, 
312, 321, 231, 213. Under S3, there are six possibilities for the isospin state of the 
three 7r's: a totally symmetric state |S), a totally antisymmetric state \A), or one 
of four mixed states |Mj) (i = 1-4). These can be defined as 

|S) = ^(|123) + |132) + |312) + |321) + |231) + |213)) , 
V6 

|Mi) ee — L(2|123) + 2|132) - |312) - |321) - |231) - |213>) , 

9 We assume that, for the indirect CP asymmetries, the CP of the final state can be fixed as for 
the decays in previous sections. Otherwise there are 2 additional theoretical parameters. 
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M- 2 ) = 


7I< |312> - 


|321) - |231> + |213» , 




M 3 ) = 


7J(-I 312 > 


- |321) + |231) + |213» , 




Ma) = 


VT2 (2|123 > 


— 2 |132> - |312) + |321) - 


|231) + |213>) 


1-4) = 


> 23 >- 


|132) + |312) - |321) + |23i; 


) - |213)) . 



(45) 

This choice of mixed states implies that two truly identical particles go in positions 
2 and 3. Under the exchange 2 -h- 3, \Mi) and |M 2 ) are symmetric, while |M 3 ) and 
|M 4 ) are antisymmetric. 

For the four B —t ixixix decays, we have: 

1. — > 7r 7r°7r°: all final-state particles are the same, which means 1 123) = 
|132) = |312) = |321) = |231) = |213). In this case, only the state \S) is 
allowed. 

2. B + ->■ 7r+7T 7r : particle 1 is tt+, particles 2 and 3 are tt°. Thus, |123) = |132), 
|312) = |213), |231) = |321). This implies that each of |M 3 ), |M 4 ), \A) is not 
allowed. 

3. B + — > 7r~7r + 7r + : particle 1 is 7r~, particles 2 and 3 are n + . Thus, 1 123) = 1 132) , 
|312) = |213), |231) = |321). This implies that each of |M 3 ), |M 4 ), \A) is not 
allowed. 

4. B% — > 7T + 7r°7r _ : we choose the order such that particle 1 is rr + , particle 2 is 7r°, 
particle 3 is 7r~. All six states are allowed. 

The amplitude for a decay with two truly identical particles has an extra factor of 
l/\/2; with three truly identical particles, the factor is 1 / a/6. 

The six elements of S3 are: I (identity), Pyi (exchanges particles 1 and 2), 
P13 (exchanges particles 1 and 3), P 2 3 (exchanges particles 2 and 3), P cy d%c (cyclic 
permutation of particle numbers, i.e. 1 — > 2, 2 — > 3, 3 — > 1), P an ticyciic (anticyclic 
permutation of particle numbers, i.e. 1 — > 3, 2 — > 1, 3 — > 2). Under the group 
transformations, \S) — > \S) and \A) — > ± \A). It is easy to see that \M\) and |M 3 ) 
transform among themselves. Writing 

\Mi)=(l) , |M 3 > = (J) , (46) 



we can represent each group element by a 2 x 2 matrix: 
I- (I °) Pi2= I J ! Pi, = 
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2 2 



1 J ' 12 ~ I 1 / ' 1 13 ~ \ 1 



1 _V3' 



-^23 ^ g _j J ) P cyclic y J/3 2 1 / ' ^anticyclic ^ ^ 1 j ' 
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Similarly, if we write 

|M 2 )^(J) , |M4> = (J) , (48) 

the 5*3 matrices take the same form, showing that \M 2 ) and |M 4 ) also transform 
among themselves. 

The above allows us to express the amplitudes for all B — > titttt decays in terms of 
diagrams. We begin with some general comments about diagrams. As an example, 
consider T\. In principle, there are six possibilities, T^ k , in which the final-state 
pions i, j, k run from top to bottom of the diagram in all permutations. Suppose 
that we want the expression for the amplitude of B — )■ 7ri7r 2 7T3 in a particular |S 3 ) 
state, and suppose that the diagram T^ k contributes to the decay. For \S 3 ) = \S), 
we define Tf: 

T f = -L (if 3 + T\ 132 + Tf 2 + T 321 + T 231 + T 213 ) . (49) 
Each T^ k leads to Tf in the amplitude. For \S 3 ) = \A), we have 




Again, each T^ k leads to in the amplitude, with a coefficient of 1 (—1) if ijk is 
in cyclic (anticyclic) order. 

For the mixed states, one has to take into account the fact that, under group 
transformations, there is |Mi)-|M 3 ) and |M2)-|Mi) mixing. In order to illustrate 
how this is done, we focus first on the M1/M3 sector. We define 

rpM\ ^ ^2^~'123 _j_ 2^~'132 rp312 rpS21 rp231 rp213^ 

T M 3 _ 1 /_ T 312_ T 321 +T 231 +T 213A ^ 

a/4 ^ ' 

Suppose I S3) = I Mi). The contribution to the amplitude of B — > 7Ti7r 2 7r3 is [M x 
(T x 1 , T x 3 ) T ] up per component, where M is the matrix representing the £3 group element 
which transforms ijk to 123 [Eq. (|47|) ]. In general, this is a combination of T x 1 and 
T^ 3 (though the T^ 3 component can be zero if M = I or P23). Factors of —1 for 
each u and l/v2 for each 7r° must also be included. If \S 3 ) = IM3), the contribution 
to the amplitude is [M x (T x Ml , Ti l3 ) T ]i ower component- This can be applied analogously 
to the M2/M4 sector, where we define 
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The entire procedure holds for all diagrams^. 

With these rules, we can now work out the amplitudes for all decays. We begin 
first with |53) = 15). The amplitudes are 



v/3 



A i B d -> T 7T 7T )| S) 



+ 



y/2A(B+ -)• 7T + 7r 7r°)| S> 



+ 



i=A(fi + ^7r-7r + 7r + )| S> 
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Cf -C 2 S +P, 

1 

3 

Pf + Cf + P* 
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1 pC„s 
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tc ' o EW1 r EW2 <^ r EWl o r EW2 
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(53) 



where P = Pi + P 2 and all amplitudes have been multiplied by e l/3 . 
For the M1/M3 sector, the amplitudes are 



V2A(B+ ^ n + n\°) lMl) 

- P5 1 + Vzp5 3 



rpM 1 V3^ M , 

2 1 2 1 2 1 -1- 2 2 2 2 
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3P t f 1 - 3\/3P t f 3 



10 Whcn applied to the decays in the previous sections, this method produces the same amplitude 
decomposition as when we used the simple rule of adding a minus sign to diagrams in which the 
identical particles are exchanged (e.g. in £> — s- Kirn or KKK). 
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cpMi , ./qpM 3 opMi i o,/qpM 3 
or EWl v ° r EW\ ° r EW2 ' ° v °- r EW/2 

p C,Mi _ , /qpC,M 3 _ pC,Mi , ./^pC,M3 
.0_- " _ " _ " 



r EWl ° v or EWl - r EW / 2 v °- r Eiy2 J ) 

2V6 A(B° d -> 7r + 7r 7r-)| M3> = [-3T x Ml + v^Tf 3 - 4>/3Tf s + SCf 1 + v^Cf 3 
+ 3C 2 Ml + v^Cf 3 + - SVSP™ 3 ] e~ ia + [-3P^ + 3^ 



- 3^^ 



pMi 
, pCMi 



/qpC,M 3 _ r pC,Mi /^ p C,M 3 
V or EWl or EW2 v or EW2 



(54) 



For the M 2 /M 4 sector, the amplitudes are 

V2A(B + -> 7T + 7T 7r )| M2> 
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(55) 



Finally, for IS3) = \A), we have 



->■ 71 71 71 J I A) 
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(56) 
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Now, the final state has isospin l<g}l<g>l = 0ffilffil©l©2©2©3. Given 
that the B- meson has I — ~ and the weak Hamiltonian has AI — ~ or |, there are 
9 paths to the final state. We therefore expect four relations among the 13 decay 
amplitudes. This is indeed what is found: 

V2A(B° d 7r°7rV)| 5> = ~V3A(B° d -> n + n\-) ]s) , 
2A(B + -> 7T + 7r 7r°)| 5> = -A{B + ->■ 7r-7r+7r+)| S> , 

|A(SS -> TrW"),^) + ^A(B° d -> TrW")^, = 

A(£?+ _> 7r + 7r vr )| A/l) - A(£?+ -> 7T-7T + 7r + )| Ml> , 

-> 7r + 7r°7r-) |M2> + ^A(B° d -> TrW"),^ = 

-)« 7r+7r 7r°)| A/2> - -)« 7T-7r + 7T + )| M2> . (57) 

These relations can also be found using the Wigner-Eckart theorem. 

In passing, we note that, within the SM, the final state with / = 3 is unreachable. 
This then provides a test of the SM. Applying the method of Ref. [IS] to B — >■ mrn, 
one can distinguish the various isospin final states. One can then look for a state 
with 1 = 3. If one is observed, this will be a smoking-gun signal of new physics. 



7.1 Dalitz Plots 

Above, we presented the amplitudes for each of the six S 3 states of B — > mm. The 
obvious question is then whether these states can be distinguished experimentally. 
Below we show that this can indeed be done. 

Consider the decay B d — > 7r + 7r°7r~. The Dalitz-plot events can be described by 
s + = (p n o +p 7r +) 2 and s_ = (p^o + p n -) 2 , so that the decay amplitude, A / l(s + ,s_ 
can be extracted. We introduce the third Mandelstam variable, so = (p-k+ + Pn- 
It is related to s + and s_ as follows: 

s + + s_ + s = m\ + 3ml . (58) 

The totally symmetric SU(3) decay amplitude is then given by 

\S) = -^[A^(s+,s_)+A4(s_,s + )+M(s + ,s ) 

+ M(s ,s+) + M(s 0l s_) + M(s_,s )} . (59) 



|Mi) = — L [2M(s+,sJ) + 2M(s_,s + ) - M(s+,s ) 



Also, 



'12 

- M(s ,s+) -M{a Q ,sJ) -M(s-,s )] . (60) 

The remaining S3 states can be found similarly. The method is similar for the other 
B — > 7T7T7r decays. 
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7.2 Weak-Phase Information 



In the previous subsection we showed how all six B — > tttyty S3 states can be ex- 
perimentally separated. It may then be possible to extract clean information about 
weak phases. (Note: by measuring the 5*3 states, one fixes the CP of the final states, 
which makes the indirect CP asymmetries well-defined.) 

Consider [S3) = \A). Here there is one decay, which yields three observables: 
the branching ratio, the direct CP asymmetry, and the indirect CP asymmetry of 
Bd 7r+7r ° 7r_ ||A}- The amplitude is expressed in terms of two effective diagrams: 
A(B® — > 7T + 7i°'n-^)\A) = Die~ ia + D 2 , which has four theoretical parameters - the 
magnitudes of -D ii2 , the relative strong phase, and a. Since the number of theoretical 
unknowns is greater than the number of observables, one cannot obtain a. Things 
are similar for \Ss) = \S). Due to the first two relations in Eq. (1571) . there are only 
two independent decays, yielding 5 observables. However, there are 8 theoretical 
parameters, so that, once again, a cannot be extracted. 

Things are different for the case of mixed states. Consider the M1/M3 sector. 
There are four decays: (1) B + -> 7r+7r 7r°|| Ml> , (2) B + ->■ 7r-7r+7r + || Ml> , (3) B° d -> 
7r + 7r°7r~|| Ml }, (4) B® — > 7r + 7r°7r~|| A f 3 }. These yield 10 observables: 4 branching ratios, 
4 direct CP asymmetries, and 2 indirect CP asymmetries (of B® — > 7r + 7r 7r^||5 3 ), 
5*3 = Mi, M3). The four decay amplitudes all have the form D\^e~ ia + D 2j i, i — 1-4. 
The Dij are related to one another by the third relation in Eq. fl57|) . as are the 
Z?2,i- The amplitudes are thus a function of 6 effective diagrams, resulting in 12 
theoretical parameters: 6 magnitudes, 5 relative strong phases, and a. Since the 
number of theoretical unknowns exceeds the number of observables, a cannot be 
extracted. However, if one assumes that the hierarchy of Eq. (|44p holds for three- 
body decays, all EWP diagrams can be neglected, to a good approximation. In this 
case, all the D 2i i are proportional to P^ 1 — V3P t A c l3 - There are thus only 4 effective 
diagrams, which yield 8 theoretical parameters. Now the number of theoretical 
unknowns is smaller than the number of observables, so that a can be obtained 
from a fit to the data. (It is not even necessary to measure all 10 observables. A 
difficult-to-obtain quantity, such as the direct CP asymmetry in B + — > 7r + 7r°7r°|| Ml ), 
can be omitted.) A similar method holds for the M2/M4 sector. The error on a can 
be reduced by comparing the two values found. 

Now, it must be conceded that the above analysis is quite theoretical - it is far 
from certain that this can be carried out experimentally [and there is an uncertain 
theoretical error due to the assumption of Eq. (I44j) ]. Still, it is interesting to see 
that, in principle, clean weak-phase information can be obtained from B — > mrn, 
or, more generally, from B — y M1M2M3 decays. 
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8 Conclusions 



In this paper, we have expressed the amplitudes for B — > M 1 M 2 M 3 decays (M, is 
a pseudoscalar meson) in terms of diagrams, concentrating on the charmless final 
states Kirir, KKK, KKn and tttttt. The diagrams are similar to those used in two- 
body decays: the color-favored and color-suppressed tree amplitudes T and C, the 
gluonic-penguin amplitudes P tc and P uc , and the color-favored and color-suppressed 
electroweak-penguin (EWP) amplitudes P E w and P% w . Here, because the final 
state has three particles, there are two types of each diagram, which we call Ti, T 2 , 
Ci, C 2 , etc. 

We have also demonstrated how to use the Dalitz plots of three-body decays to 
separate the decay amplitudes into pieces which are symmetric or antisymmetric 
under the exchange of two of the final-state particles. This is useful for any decay 
whose final state contains identical particles under isospin. If the relative angular 
momentum of the two particles is even (odd), the isospin state must be symmetric 
(antisymmetric). These two possibilities can be distinguished experimentally. 

The main advantage of a diagrammatic analysis is that the approximate rela- 
tive sizes of the diagrams can be estimated. For example, there are annihilation- 
and exchange-type diagrams which contribute to these decays. However, these are 
expected to be negligible, and are not included in our analysis. Previous studies of 
three-body decays were carried out using isospin amplitudes, and gave exact results 
for the symmetric or antisymmetric states. On the other hand, the (justified) neglect 
of annihilation-type diagrams can modify these results, and can lead to interesting 
new effects. 

As an example, consider B — > KKK, which consists of four decays. For the case 
where the two K's are in a symmetric isospin state, the Wigner-Eckart theorem 
gives a single relation among the four amplitudes. However, when the amplitudes 
are written in terms of the non-negligible diagrams, it is found that this relation 
actually consists of two equalities, and this leads to new predictions of the standard 
model (SM). Present data allow us to test one of these equalities, and we find 
agreement with the SM. In the same vein, B — > KKK decays can be written in 
terms of five isospin amplitudes. The diagrammatic analysis shows that, in fact, 
only four of these are independent - two of the isospin amplitudes are proportional 
to one another. 

Another consequence of the diagrammatic analysis has to do with weak phases. 
The CP of a three-particle final state is not fixed, because the relative angular mo- 
menta are unknown (i.e. they can be even or odd). For this reason, in the past it 
was thought that it is not possible to cleanly extract weak-phase information from 
three-body B decays. In this paper, we demonstrate that this is not true. Using 
the diagrams, we show that it is possible to cleanly measure the weak phases in 
some decays, given that it is experimentally possible to distinguish different symme- 
try combinations of the final-state particles. We explicitly give methods for KKn 
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and tctttc, and note that the the procedure for Kirir is presented separately. Ways 
of cleanly extracting the CP phases from other three-body decays will surely be 
suggested. 

There are thus a number of interesting measurements that can be carried out 
with B — > M1M2M3. LHCb is running at present, and the super-5 factories will run 
in the future. Hopefully, these machines will provide interesting data on three-body 
B decays. 
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